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Abstract 

Hydrodynamic forcing plays an important role in shaping the dy- 
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namics of marine organisms, in particular of plankton. In this work 
we study the planktonic biological activity in the wake of an island 
which is close to an upwelling region. Our research is based on nu- 
merical analysis of a kinematic flow mimicking the hydrodynamics in 
the wake, coupled to a three component plankton model. 

We use parameter values of relevance for the Canary wake, and the 
main results for a realistic range of parameters in this area area are: 
a) Primary production is enhanced in the region of the wake opposite 
to the upwelling zone, b) There is a strong dependence of the pro- 
ductivity on the inflow conditions of biological material entering the 
wake transported by the main current. Finally c) we show that under 
certain conditions the interplay between wake structures and biolog- 
ical growth leads to plankton blooms inside mesoscale hydrodynamic 
vortices that act as incubators of primary production. 

Keywords: plankton; island wake; primary production; upwelling; 
vortex dynamics 
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1 Introduction 



Understanding the influence of hydrodynamic motions on the growth, pro- 
ductivity and distribution of marine organisms, especially in the context of 
plankton dyna mics, is a major challenge recently addressed from a variety 



of perspectives (|Mann &: Lazier!. 



1991 



Penman fc Gargett 



1995 



Abraham 



1998 



2002 



Peters &: Marrase 



Martin 



2000 



Karolyi et al 



2000 



Lopez et al. 



2001a 



Martin et al. 



20031 ). Vertical transport processes of nutrients are recognize d 



19951 ). 



as key factors controlling plankton productivity (iDenman fc Gargett 
In particular, upwelling areas in the world's oceans are of fundamental im- 
portance for the growth of phytoplankton which is the base of oceanic food 
webs. They are characterized by nutrient rich waters coming to the sur- 
face from depths of over 50 meters. Nutrient enrichment enhances phyto- 
plankton growth close to the upwelling regions, giving rise to an increase 
in zooplankton and fish populations in the area. More rec ently, the impor 



tance of hor i zontal fluid motion has also been pointed out ([Abraham!. 



Lopez et al.. 


2001b 


2003; 


Martin 


. 2003) 



Hernandez- Garcia et al.l . [2002 



1998; 



Hernandez- Garcia et al.l . 



20031 ). Mesoscale stirring redistributes and mixes 



nutri ents laterally, giving rise also to enhanced pr oductivity (jMartin et al. 



plankton anc 



20021 ). or to bloom initiat ion (IReigada et al. 



petition and coexistence (jKarolyi et al. 



2000 



2003) , and affects sp ecies com 



Bracco et al. 



20001 ). Satellite 



images illustrate the interaction between horizontal mesoscale motions and 
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plankton dynamics. 

Vertical upwelling and strong mesoscale activity occur simultaneously in 
several places of the globe. A stronger impact and a high complexity of the 
physical-biological interactions are expected there. Some of these areas are 
the Benguela zone, the Humboldt Current, or the Canary islands. 

Though the phenomena we discuss are rather general, we illustrate them 
by using the Canary islands, which are close to the northwestern African 
coast, as a specific example. There, upwelling occurs at the African coast 
because of Ekman pumping induced by the dominant winds, and in addition, 
the Canary islands constitute an obstacle for the main ocean current in the 
area, flowing from Northeast to Southwest, originating a strong mesoscale 
hydrodynamic activity in their wake. The interaction between the vortices 
in the wake and the Ekman flow transporting nutrient-rich waters from the 
coastal upwelling seems to be at the heart of the observed enhancement of 
biological production in the open Atlantic ocean close to the Canary region. 
Motivated by this situation, the aim of this paper is to study, in a more gen- 
eral framework, the role of wake vorticity in redistributing upwelled nutrients 
and influencing phytoplankton growth. 



To this end we combine the kinematic model flow introduced in 



Sandulescu et al. 



( 120061 ) with a simple model of a Nutrient-Phytoplankton-Zooplankton (NPZ) 



trophic chain, and study the impact of the flow characteristics on the bio- 
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logical dynamics, particularly on the primary production (PP). We will use 
mainly parameter values of relevance in modeling the Canary wake, but we 
expect our results to have broader application. Only horizontal transport is 
explicitly taken into account in the flow, the upwelling is modelled as a source 
term in the nutrient equation. We address questions such as (i) whether the 
island wake is a barrier for the upwelled nutrients, or (ii) if rather the gener- 
ated stirring mixes nutrients into poorer waters so that primary production 
is enhanced, or (iii) what is the impact of the presence of vortices and other 
wake structures on biological activity. Our main results are, on the one hand, 
that for a range of parameters which is realistic in the Canary area, primary 
production is enhanced in the part of the wake opposite to the upwelling 
zone. That is, the wake is not a barrier confining the region of high nutrients 
and plankton growth. On the other hand, there is a strong dependence of 
the productivity, and of the role of the vortices, on the inflow of biological 
components entering the wake due to transport by the main current. In some 
situations the vortices in the wake act as an incubator whose sole presence 
is enough to greatly enhance biological productivity in poor waters entering 
the region. 



Recently, in 



Sandulescu et al. 



(120071 ) we have studied the same model but 



attending to the role of the different time scales present in the system. We 
find there that the long residence times of nutrients and plankton close to the 
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island, and the confinement of plankton inside vortices are important factors 
for the appearence of localized plankton blooms. These studies complement 
the final part of this paper where we focus on the role of vortices on biological 
production. 

In the next section we present our general modeling framework, present- 
ing the velocity field, the plankton model, and the boundary (environmental) 
conditions. In section [3] we present our results, organized in two subsections 
that contain our studies of concentrations of nutrients, plankton, and pri- 
mary production in relation to hydrodynamic and inflow conditions, and 
another one in which we analyze the plankton content of vortices. Section H] 
summarizes our conclusions. 

2 Modeling framework 

Figure [I] shows our two-dimensional model domain. The main current flows 
from left to right, passing by the circular obstacle, which models the pres- 
ence of an island, and giving rise to mesoscale vortices in its wake. Vertical 
hydrodynamic motion is not explicitly considered, but its effect on nutrient 
upwelling is modelled by a source of nutrients (shown as a small box in the 
upper zone of the domain). The associated Ekman flow points towards the 
interior of the domain. Our focus of study will be the lower part of the wake, 
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the region A s , marked with a box in the lower part of figure [TJ A discussion 
about the ch oice of a different A s and its negligible effect in our results can be 



consulted in 



Sandulescu et al 



(120061 ). Our primary objective is to determine 
if nutrient input from the upwelling region, which is on the opposite part of 
the wake, may enhance the biological activity in this region, and to elucidate 
the role, as barriers or as transporters, of the wake and of the vortices present 
in it. In addressing this goal, we realize the importance of the contents of the 
water transported towards our domain by the main flow. In the context of 
the Canary islands situation, that will guide our selection of parameter val- 
ues, the upper part of figure [T] represents the African coast, with the coastal 
upwelling. The obstacle is the Canary archipelago, more particularly the 
Gran Canaria island which plays an important role for the emergence of the 
vortices in the area. The main current is the Canary current, flowing from 
northeast to southwest. In this context our analysis may be of relevance to 
discuss enrichment of the open ocean beyond the Canary wake by input of 
coastal waters. More generally, it illustrates the interplay between transport, 
stirring, and biological dynamics. 



2.1 The velocity field 

We briefly introduce the velocity field used in this study. It is essentially the 
horizontal incompressible flow, derived from a time periodic streamfunction 
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Figure 1: The computational domain with a snapshot of the velocity field. 
Spatial coordinates, x and y, are in units of the island radius. 
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of period T c , proposed by 



Jung et al. 



(119931 ) to model kinematically the vor- 



tex street behind a cylinder at moderate Reynolds numbers, but modified to 
include a velocity component pointing towards the domain interior that mim - 



ics the Ekman flow associated with the upwelling (Sandulcs cu et al. 



Sandulescu et al. 



2006|) 



2006 . We next 



Its technical description can be found in 
describe it qualitatively. 

There are a maximum of two vortices simultaneously in the system. They 
are of opposite vorticity sign but their maximal vortex strength denoted by 
w is equal. They are created behind the circular obstacle with a phase differ- 
ence of half a period, T c /2. Each of the vortices travels a distance along the 
x direction for a time T c and finally disappears. Then the process repeats 
periodically again. Since real oceanic flows are never perfectly periodic we 
add some randomness to the vortex trajectories. Instead of moving along 
straight horizontal lines, they are subjected to a stochastic transversal Brow- 



Sandulescu et al 



toom ). 



nian motion of small amplitude (see details in 

The main background flow moves in the positive horizontal direction with 
a speed uo, and the Ekman drift, which is intended to model the flow from the 
coast towards the ocean interior, is introduced by considering an additional 
velocity ue in the y direction acting in the region with x coordinate larger 
than the island radius r, i.e. just behind the island (x > r) (this additional 
component ue is indicated in[fl but not added to the velocity field plotted 
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there). The circular obstacle, which is considered as a model island, has a 
radius r. 

To adapt this general setup of the velocity field to a realistic and more 
specific situation we c hoose parameter values which are guided by the values 



in the Canary wake (ISandulescu et al. 



20061 ): r = 25 km, u = 0.18 m/s, 



ue = 0.02 m/s, and T c = 30 days. We consider two situations for the vortex 



strength in the wake. Previous results (ISandulescu et al. 



20061 ) indicate that 



the wake entrains water from one side of the island towards the other in form 
of filaments for high values of the vortex strength which are realistic in the 
Canary area (w = wh ~ 55 x 10 3 m 2 /s), but that it acts as a barrier to 
transport when the vortices are weak (say w = wl ~ w#/20 = 2.75 x 
10 3 m 2 /s). We will analyze the plankton dynamics under these two vortex 
strengths, Wl and Wh, and also at intermediate ones. Despite the smallness 
of ue, it is larger than the minimum neede d to observe a transitio n from no 



transport to transport when increasing w (ISandulescu et al. 



200J). 



2.2 The NPZ model 



Our descri ption of the plankton 



popu lation dynamics is based on a model de- 



veloped by lOschlies fc Garconl (jl999l ). It describes the interaction of a three 
level trophic chain in the mixed layer of the ocean, consisting of nutrients N, 
phytoplankton P and zooplankton Z, whose concentrations evolve in time 



10 



with the following NPZ dynamics: 



dN 
~~dt 



F N = ^ N -(3 



k N + N 
arjP 2 r 



P 



dP 
~dt 
dZ 
~dt 




a + 7]P 2 




Z + fi P P + fi Z Z 2 



Z ^ a + fjP 2 



arjP 2 



(1) 



The dynamics of the nutrients includes three different processes. There 
is a nutrient supply given by $jv = S(x,y)(N — N) due to vertical mix- 
ing. S gives the inverse of the time scale for the nutrients to relax to the 
nutrient concentration iVo below the mixed layer. Therefore S is the pa- 
rameter accounting for the vertical nutrient supply due to upwelling. We 
take S(x,y) = Si = 0.00648 day^ 1 outside of the upwelling region and 
S(x,y) = Sh = lOOSi = 0.648 day -1 in the nutrient-rich upwelling area 
identified in figure [TJ The nutrients are consumed by the phytoplankton 
according to a Holling type II functional response. The last three terms in- 
side the parenthesis of the nutrient equation denote the recycling of a part 
of all dead organic matter. The phytoplankton grows upon the consump- 
tion of the nutrients, but its concentration is decreased due to grazing by 
zooplankton and to natural mortality. The grazing enters as a growth term 
for the zooplankton concentration with an efficiency factor 7. Zooplankton 
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parameter 


value 




0.66 day -1 


V 


1.0 (mmol N m -3 )~ 2 day -1 


7 


0.75 


a 


2.0 day" 1 


s t 


0.00648 day -1 (nutrient poor) 


s h 


0.648 day -1 (nutrient rich) 


k N 


0.5 mmol N m -3 


H N 


0.2 


fi P 


0.03 day - 1 


Hz 


0.2 (mmol N uT 3 ) -1 day -1 


N 


8.0 mmol N m -3 



Table 1: List of parameters used in the biological model 



mortality is assumed to b e qua dratic. Additional det ails can be consulted in 



Oschlies fc Garconl ( 



1999) and 



are taken from iPasquero et al 



Pasauero et al. 



(120051 ). The parameters used 



(120041 ) and presented in the Table [TJ 



The primary production, the rate at which new organic matter is pro- 
duced, is given by the growth term in the phytoplankton dynamics 



PP = (3 



N 



k N + N 



-P. 



(2) 



The dynamics of this food chain model is studied in detail in 
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Edwards & Brindley 



(119961 ) and iPasquero et al.l (120041 ). Depending on the parameters of the 
model, it exhibits stationary or oscillatory behavior in the long-term limit. 
The chosen parameter values lead to a steady state. Using the values from Ta- 
bleland fixing the vertical mixing to the lower value S = Si = 0.00648 day -1 
we obtain a steady state that will be called the ambient state: N am b = 0.185, 
Pamb = 0.355 and Z am b = 0.444 mmol N m~ 3 . Thus, in the nutrient poor 
region occupying most of the domain the ambient primary production in 
steady state is PP am b = 0.0633 mmol N m~ 3 day -1 . In the upwelling region, 
S = Sh = 0.648 day -1 and the steady state that would be reached un- 
der this nutrient input would be N up = 7.539, P up = 0.603 and Z up = 1.154 
mmol N m -3 , and the primary production associated with these values would 
be PP up = 0.373 mmol N m -3 , i.e. nearly 6 times ?? am t,. 



2.3 Complete model and input conditions 

The coupling of the biological and the hydrodynamic model yields an advection- 
reaction system. We add also an eddy diffusion process acting on plankton 
and nutrients concentrations with diffusion coefficient D to incorporate the 
small scale turbulence, which is not expli citly taken in to account by the 



Okubol (ll97ll ) prescriptions, we 



large scale velocity field used. Following 
take D pa 10m 2 /s, corresponding to spatial scales of about 10 km at which 
flow details begin to be absent from our large scale flow model. Thus our 
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complete model is given by the partial differential equations: 




F N + DV 2 N, 



F P + DV 2 P, 



F z + DV 2 Z, 



(3) 



with the biological interactions Fn, Fp, and Fz from Eq. (JTJ), and the velocity 
field v(x,y,t) described in subsection 12.11 This system is numerically solved 
by means of a semi-Lagrangian algorithm on a grid. Additional details of 
the integration algorithm are reported in the Appendix |A] 

Since we are studying an open flow, inflow conditions into the left part of 
the domain should be specified. It turns out that the influence of inflow con- 
centrations is rather important and we present here two cases that exemplify 
the two main behaviors we have identified: In the first one fluid parcels enter 
the computational domain with the ambient concentrations N am b, P am b, and 
Zamb- This corresponds to the steady state for S = S[, and represents the 
situation in which the exterior of the computational domain has the same 
properties as the part of the domain without upwelling. This input condition 
will allow us to focus on the interaction between the upwelling water and the 
main part of the domain containing the wake. In our second situation fluid 
particles transported by the main flow enter the domain from the left with a 
close to vanishing content of nutrients and plankton, corresponding to a bio- 
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logically very poor open ocean outside the considered domain. To be specific, 
we take N L = 0.01N amb , P L = 0.01P amfe and Z L = 0.01Z amb . Primary pro- 
duction in the inflow water is very low: PPl = 8.6 x 1CT 6 mmol N m~ 3 day -1 . 
This is more than 7000 times smaller than PP amb . Since those concentrations 
are very low, we take into account that fluctuations may be important by 
adding to each of the concentrations (N L , P L , Z L ) of each fluid parcel enter- 
ing the system an independent random amount of about ~ 5% of the inflow 
concentration. In this second situation there is mixing between three types 
of water: the 'ambient', the 'upwelled', and the 'inflow' ones. It turns out 
that the interaction between inflow and wake will be the responsible for the 
interesting behavior described below. 



3 Results 

In this section we first describe the outcome of the different scenarios con- 
sidered in terms of primary production and plankton distributions, and then 
address in more detail the relation between vortex structures and plankton 
patches. The brackets <> shall denote spatiotemporal averages. 
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5 10 15 

t/T 

c 

Figure 2: The time evolution of the ratio between the primary production, 
PP, spatially averaged in A s , and the ambient one PP amb for the two inflow 
cases and two values of the vortex strength. The two upper lines are for 
ambient input concentrations, and the two lower for low inflow. Dashed- 
dotted lines w = wl, solid lines w = wh- 
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3.1 Primary production and plankton dynamics 



Sandulescu et al. 



fl2006h is 



We stress that one of the main observations in 
that, for the flow parameters used here, there is a qualitative change in the 
transport behavior at vortex strength w c ~ u;#/10: For weaker vortices, a 
plume of passive tracers released from the location of our upwelling area de- 
velops in the direction of the main flow with a slight transverse displacement 
due to the Ekman flow ue but remaining far from our study region A s . The 
wake acts here as a barrier to transport. For w > w c , however, the plume 
becomes a filament that is entrained by the vortices, so that it crosses the 
wake and reaches A s . Note that wl < w c < wh- 

The observed behavior of our biological model when ambient concentra- 
tions are used at the inflow reflects directly this transport behavior: The two 
upper lines in figure [2] show the time evolution of the productivity averaged 
over the region A s . The dashed-dotted line with a nearly constant value 
PP/ PP amb « 1 is obtained for w = Wl < w c . The upwelling plume fertilizes 
the upper part of the computational domain, where higher concentrations of 
plankton are observed, but the lower part of the wake is unaffected by this 
and keeps its low ambient productivity value PP am b nearly constant. When 
w = wh > w c (upper solid line) productivity becomes enhanced with respect 
to its ambient value. It undergoes roughly periodic oscillations reflecting 
the periodic motion of the nutrient filament entrained by the vortices. The 
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central column in figure [3] displays the phytoplankton spatial distribution 
at different time instants. A filament of high phytoplankton concentration 
appears in the system, sitting basically on top of the high nutrient filament 
(not shown) emerging from the upwelling and being entrained by the vor- 
tices. Zooplankton and primary production are also distributed in the same 
way. 

Figure H] shows the average primary production (averaged over A s and 
then averaged in time) as a function of the vortex strength w in the range 
[0.025iOff, Wh\- As anticipated, a transition from essentially no enrichment 
by the upwelling to an increasing primary production occurs around w = 
w c ~ O.Iwh, confirming a direct influence of the physical transport process 
on the biological dynamics. 

The dynamics in the low concentration inflow case is very different. For all 
values of w considered, the average primary production in A s is smaller than 
the ambient one. This can be understood from the fact that fluid elements 
enter the domain with very low nutrient and plankton concentrations. 

To better understand this situation we have performed numerical analysis 
of the NPZ dynamics without hydrodynamic terms. Figure [5] shows the time 
needed for P and PP to reach their maximum values as a function of the 
initial conditions. We define / as the fraction of the ambient concentrations 
that are present in the inflow: (N,P,Z) — fx (N amb , P amb , Z amb ). With a 
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Figure 3: The Okubo- Weiss parameter (left column), the concentration of 
phytoplankton for ambient inflow (middle), and concentration of phytoplank- 
ton for the low inflow situation (right). Phytoplankton concentration is ex- 
pressed in units of P a mb, and w = wh- Snapshots taken during one flow 

period, at t/T c = 14, 14.25, 14.5, 14.75, from top to bottom. 
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mean flow of speed u = 0.18m/ s, fluid elements spend only 16 days (0.53 in 
units of T c ) inside the domain of horizontal extension lOr = 250km in figure 
[TJ We see in Fig. [5] that for / < 0.4 the maximum in P occurs later, so that 
we cannot expect considerable growth in A s outside the vortices for the value 
of / = 0.01 used in this work. But this observation is also puzzling, since 
the primary production reported in figure [2] for the low inflow case is not as 
small as the above argument would indicate: It is reduced just between 40% 
and 85% with respect to the ambient values. 

Figure [3] (right column) clarifies the mechanisms involved. The spatial 
plankton distribution is rather different from the ambient inflow case. It 
is clearly related to the vortices and associated structures. The plankton 
concentration is very low outside these objects. As in the case of ambient 
inflow, a plume with high nutrient concentration is present in the system 
due to upwelling and has a shape similar to the one in the ambient inflow 
case, which resembles the phytoplankton distribution of the central column of 
fig. [31 but here it seems to have no effect in inducing phytoplankton growth. 
The time scale for plankton growth starting from small values is much larger 
than the travel time through the computational area, thus these effects are 
observable only further downstream. Therefore, in the study area displayed 
in figures [1] and [3] the influence of the upwelling nutrient filament is masked 
by a more prominent mechanism, described below. In fact, in this low- 
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inflow case, we note that the phenomenology observed in the area remains 
qualitatively unchanged if the upwelling is removed, although quantitative 
changes occur. This indicates that the dominant mechanism in the low inflow 
situation is not the mixing of upwelling and ambient waters, as in the ambient 
inflow situation, but the interaction of the inflow with the wake. 

Figure [3] shows that phytoplankton growth in the vortices occurs after 
phytoplankton is transported into their interior by filaments emerging from 
the boundary of the circular obstacle. This complex structure - boundary 



is wel 



1993 



known from dynamical systems studies o 



Ziemniak et al. 



1994 



Karolyi et al. 



2000 



Tel et al. 



2005 



(June et al. 




Sandulescu et al. 



20061 ). and is related to the so-called unstable manifold of the chaotic sad- 
dle, the main dynamical structure in the wakes occurring in time-dependent 
two-dimensional open flows. Loosely speaking, it is the location of the fluid 
elements that take a long time to leave the proximity of the island, because 
of the complex recirculation emerging just behind the obstacle as well as the 
reduced velocities occurring near its boundary. A detailed analysis of these 
st ructures and their implic ations for the residence times has been published 



in 



Sandulescu et al. 



(120071 ). Particle release calculations allow us to realize 
that, although most of the incoming particles follow the mean flow and leave 
the system in the 16 days lapse estimated before, a fraction of them are cap- 
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Figure 4: The spatiotemporal average, < PP >, of the primary production 
in A s in terms of its ambient value, PP am b, as a function of the vortex 
strength w normalized by Wh- We plot the two inflow cases: ambient inflow 
concentrations (•) and low inflow concentrations (o). 
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Figure 5: Time needed, in units of T c , to reach the maximum value of P or 
of PP, vs /, for the NPZ dynamics without flow. As indicated in the legend, 
squares are for P and circles for PP. 
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Figure 6: The spatiotemporal average of the primary production in A s in 
terms of its ambient value < PP > /PP am b as a function of the inflow 
concentrations at vortex strength w = Wh (solid line) and w = Wl (dashed 
line). / is the fraction of the ambient concentrations in the fluid entering the 
system. 
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tured by the wake structures with residence times of about 50 to 20 days 
for wh and Wl, respectively These long residence times allow the plank- 
ton concentration to build up in the filaments emerging from the obstacle, 
which gives rise to a plankton bloom later downstream, when the filaments 
are stretched and rolled up by the vortices. 

Thus, recirculating structures in the island wake act as incubators that 
make fluid elements more productive before releasing them into the main cur- 
rent. It turns out that the peak values of the phytoplankton bloom in this 
low inflow case are larger than the ones under ambient inflow. This somehow 
paradoxical observation is explained by the fact that zooplankton values are 
relatively high under equilibrium ambient conditions, so that grazing con- 
trol of the phytoplankton population is rather effective. By contrast, in the 
low inflow case zooplankton and thus grazing control is essentially absent. 
Zooplankton concentration begins to build up only when phytoplankton con- 
centration has already reached larger values. Therefore it is responsible for 
the end of the bloom further downstream, but high phytoplankton values are 
attained before that. 

Figure [2] (two lower curves) shows the time evolution of the primary 
production under the low inflow conditions. Even for w = w^, for which the 
wake acts as a barrier blocking nutrient fertilization of A s from the upwelling, 
primary production shows an oscillating behavior, reflecting the oscillations 
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of the wake structure which is the responsible for the plankton growth. Figure 
H] shows the increase in primary production by increasing the vortex strength 
w starting at a value of w = w c ~ O.lw. In the range of w considered there 
is an increase in primary production by a factor of about 2.17, larger than 
the factor 1.2 of increase attained under ambient inflow. 

So far we have described two very distinct inflow situations and studied 
the impact of vortices by varying the vortex strength. We now fix the vortex 
strength w to the high, wh, and low, wl, values, and describe the primary 
production behavior for intermediate inflow cases in figure [6] (solid line is for 
Wh and dashed line for wl). The inflow concentrations are now varied in 
terms of / which, as already mentioned, gives the fraction of the ambient 
concentrations that are present in the inflow. For / = 0.01 we have the low 
inflow conditions considered before, and / = 1 corresponds to ambient input. 
We see, for both values of w, an increase of the average production in A s with 
the biological content of the inflow, which is rather fast until / « 0.2. 

Maximum values of the averages of P and PP in A s are shown in Fig 
for different initial conditions, with and without flow (i.e. in the spatially 
homogeneous case). It is seen a contrasting behavior remarking the fact that 
the flow is not simply redistributing biological material, but that it modifies 
also the spatial averages and then the total amounts of substances, as well as 
the qualitative dependence on /. The large overshoots of the homogeneous 
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case for low / only occur inside vortices in the presence of flow, as will be 
discussed in the next section, so that the spatially averaged values of P and 
PP are lower (see also Fig. [2]) . When approaching / = 1 (the ambient inflow 
case), the relative increase in the presence of flow arises by advection from 
the upwelling plume. 

For completeness, we present in FigJH] the effect of reducing the effective 
diffusion D. Physically this would mean that we decrease the intensity of 
small scale turbulence. We see that in general primary production is slightly 
reduced. The increase which is observed at low values of w in the low inflow 
case is due to the lack of dilution of the filaments which emerge from the 
cylinder boundary. For the ambient inflow situation a small decrease of 
productivity for smaller D is seen only at large vorticity, because distributions 
in A s are rather homogeneous otherwise. 



3.2 Vortices and plankton distribution 

It is well known that vortices are responsible for a large part of the trans- 



1998: 


Peleeri et al.. 


2005; 


Martin et al.. 


2002 



Barton et al. 



2002). They influence biological dy- 



namics, and most of the studies have focused on the effect of the relatively 
large vertical motions induced by their cores. Here we focus instead on hori- 
zontal processes. In this section we consider the case w = Wh, where strong 
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Figure 7: Maxima of the average in A s of P and PP (normalized with the 
corresponding ambient values) vs /. The circles are for the P concentrations 
and the squares for PP. Data points joined with solid line are computed 
without flow, and the ones joined by dashed line are computed with flow. 
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Figure 8: Average value of PP in A 5 computed for a small value of the 
diffusion coefficient, D s = 0.1D, over the one obtained for D, versus w/wh- 
Filled circles correspond to the low inflow case, while squares to the ambient 
inflow one. 
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vortices are present in the system, and characterize the plankton distribu- 
tions relative to vortex positions for the two different inflow conditions, that 
highlight the two different primary production enhancement mechanisms dis- 



(Okubo. 


f970; 


Weiss. 


1991) 



W (a precise definition is included in the Appendix [B]) to identify in an 
objective way the interior and the exterior part of vortices. Flow regions 
with W < are vorticity dominated, and can be identified as the inner 
part of vortices. Regions with W > are strain dominated and outside 
vortices. The leftmost column of figure [3] displays the W values, showing 
clearly the position of the vortices. These positions can be correlated with 
the phytoplankton distributions displayed in the other columns. 

Figure M shows the spatial average of phytoplankton concentrations and 
primary production inside each of the two vortices and outside them. These 
regions are identified with the help of the Okubo- Weiss parameter W. The 
primary production time series is qualitatively similar to the phytoplankton 
one, although slightly shifted towards earlier times. This is so because PP 
contains the influence of the nutrient dynamics, whose temporal evolution 
anticipates the phytoplankton one. The zooplankton time series (not shown) 
are also qualitatively similar but shifted towards later times. 

In the case of ambient inflow the interior of the vortices (dashed and con- 
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Figure 9: Upper panels: time evolution of the spatial average of the phy- 
toplankton concentration P (normalized with P am b) inside the two vortices: 
y < (lower vortex, solid line) and y > (upper vortex, dashed line), and 
out (outside, dot line). Left is for ambient inflow concentrations, and right 
for low inflow concentrations. Vortex strength is always w = Wh- Lower 
panels: primary production in the same locations and situations. 
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tinuous lines) contains the same quantity of plankton as the inflow, namely 
the ambient concentration. Only when additional nutrients from the up- 
welling zone are entrained we observe bursts localized in time. Thus most 
of the biological activity is in the outside area (which includes the upwelling 
zone) . This quantifies what is seen in the middle column of figure [3j plankton 
appears mainly in filaments that wind around the vortex periphery basically 
without entering them. The asymmetry observed between the content of the 
two vortices arises from the fact that, due to the different sense of rotation of 
the vortices more nutrients are transported towards the vicinity of the lower 
vortex than to the upper one (cf. figure [3]). 

The situation is rather different in the low inflow case. The range of the 
concentration oscillations is now larger, and the content of the two vortices 
oscillates in antiphase. The largest concentration values occur now inside 
vortex cores, leading to peak bloom values larger than for the ambient inflow 
situation. Minima are also smaller so that averages in regions such as A s 
give an overall smaller plankton content and primary production. Filamen- 
tal structures close to the boundary of the island transport concentrations 
towards the vortices where the species are trapped and transported down- 
stream. During this motion their concentrations are homogenized by small 
scale turbulence (modeled by the diffusion term in Eq. [3]) and the classi- 
cal dynamics that the system of equations ([!]) exhibits in a homogeneous 
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situation occurs: nutrient consumption by phytoplankton induces a large 
phytoplankton bloom which is stopped by the grazing by zooplankton, that 
also experiences growth, until all three components approach the final equi- 
librium value (N amb , P amb , Z amb ). This steady state for the vortex content 
occurs only further downstream. 



4 Conclusions 



We have presented numerical results on the biological dynamics in the wake of 
an island close to a coastal upwelling area. Parameter values were appropriate 
for the Canary Islands region but we expect our results to be of greater 
generality. 

Two different scenarios have been identified and discussed. In the first 
one, occurring when the region outside the focus area has properties similar 
to it, we have identified an enrichment mechanism of one side of the wake 
by nutrients upwelled on the other side. It occurs at sufficiently high vortex 
strength of the vortices present in the wake. Vortices entrain water from one 
side of the island in the form of filaments that are transported across the 



wake. Filaments o: 



area (IBarton et al 



this type are observed in satellite images of the Canary 



1998 



Pelegri et al. 



20051 ). When the vortex strength is 



low, the wake acts rather as a barrier that blocks transport from one side 



33 



to the other. This scenario is a direct translation o f the behavior of passive 



tracers under similar flow (jSandulescu et al 



20061 ). We have also observed 



that a large decreasing of the eddy diffusion coefficient value has no relevant 
role in the primary productivity, though in general this is slightly reduced. 

The second scenario becomes evident when the waters surrounding the 
study area are biologically much poorer. Now fertilization by the upwelling 
is not relevant, but we have identified a mechanism for primary production 
increase in the wake: The large residence times of some of the fluid particles 
in particular structures of the island wake allow them to become enriched by 
the ambient nutrient sources. Filaments from the wake structures feed this 
enriched water into the vortices, and the nonequilibrium plankton dynam- 
ics there leads to strong plankton blooms confined inside the vortices. The 
biological significance of hydrod ynamical structures in the wake of obsta- 



cles has been re cognized before (jKarolyi et al. 



2000 



Tel et al 



Scheming et al 



2000 ; 



20051 ). but in these cases the relevance was associated with their 
complex geometric structure that allowed fine intertwining of filaments con- 
taining different species or substances. The mechanism presented here seems 
to be different and simply associated with large residence times in the wake, 
leading to a kind of incubatory effect. A carefull study of th e role of the differ- 



ent tim escales present in the system has been presented in 



Sandulescu et al 



(120071 ). We expect this mechanism to be at work in many types of island 



34 



wakes, even if they are not associated with upwelling systems. 
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A Numerical Algorithm 

The coupling of the biological and the hydrodynamic model yields an advection- 
reaction-diffusion system. This system is numerically solved by a semi- 
Lagrangian algorithm. The concentrations of N, P and Z are represented 
on a grid of [500 x 300] points. The three processes, advection, reaction and 
diffusion, are performed sequentially as follows: 

• Advection: Every point of the grid is integrated backwards in time 
with the velocity field for a time step dt. In this way we obtain the 
past position of the chosen point, which typically is not on the grid. 
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• Reaction: the concentration of N, P, or Z at the position obtained 
in the advection step is computed using a bilinear interpolation of its 
corresponding nearest points on the grid. Then these concentrations are 
the starting values for the integration of the reaction dynamics forward- 
in-time for a time step dt. In this way, we obtain the concentration of 
N, P and Z in the original grid point. 

• Diffusion: Once we have the values of the concentrations on the grid, 
the diffusion step is performed following an Eulerian scheme. However, 
the interpolation in the reaction part induces a numerical diffusion of 
the order D n oc dx 2 /dt. It is therefore desirable to e nsure th at the real 



197l|)£> 



diffusion, which is given by the Okubo estimation (lOkubo . 
10m 2 /s, is larger than this numerical diffusion. Moreover, the stability 
condition for the Eulerian diffusion step requires that Ddtd/dx 2 < 1, 
where dtj denotes the time step for the diffusion part. These two 
inequalities are fullfilled if the diffusion time step dtd is smaller than dt. 
According to this condition we chose dtd = dt/ '10 in our algorithm. This 
implies that the algorithm makes ten steps of diffusion dt^ after every 
step of advection and reaction. Expressed in units of T c = 30 days for 
time and r = 25 km for space, the dimensionless numerical values of the 
parameters are dx = 0.02, dt = 0.01, dt d = 0.001, and D = 0.041472. 
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B The Okubo- Weiss parameter 



The Okubo- Weiss parameter W 



Okubol f|l970h 



Weiss! (jl99ll ) is a quantity 



used to distinguish areas in which the flow is dominated by vorticity from 
those areas where the flow is strain dominated. It is given by 



W = sl + sl-u\ (4) 

where s n , s s are the normal and the shear components of strain, and uj is 
the relative vorticity of the flow defined as: 

du dv dv du dv du 

n dx dy' s dx dy ' dx dy 

Where u and v are the horizontal components of the velocity field. We 
chose the critical threshold value to be W c = 0. For areas where W is below 
W c the flow is vortex dominated, otherwise we consider the flow to be strain 
dominated. 
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